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ABSTRACT 

New measurements of the stellar-mass deficits at the centers of luminous elliptical galaxies 
are presented. These were derived considering the following observational facts. Firstly, “core” 
galaxies, which are thought to have had their inner region depleted from the coalesence of su¬ 
permassive black holes, show an abrupt downward deviation of their inner light-profile relative 
to their outer Sersic profile. Second, “power-law” galaxies, having undisturbed prohles and no 
partially depleted core, have inner light-profiles that display no departure from the inward ex¬ 
trapolation of their outer Sersic prohle. The central stellar deficits have therefore been derived 
from the difference in flux between the ffST-observed galaxy light-profiles and the inward ex¬ 
trapolation of each galaxy’s outer Sersic profile. This approach gives flux dehcits ^0.1% of the 
total galaxy light, and mass deficits that are ^2 times each galaxy’s central supermassive black 
hole mass. These results are in agreement with the theoretical expectations of mass ejection from 
binary black hole mergers and also with popular ACDM models of hierarchical galaxy formation. 
It is also explained why this result is some 10 times smaller than current observational estimates 
of the central mass deficit, and therefore implies a merger history for giant elliptical galaxies that 
is one order of magnitude less violent than previously suggested. 

Subject headings: black hole physics — galaxies: fundamental parameters — galaxies: elliptical and 
lenticular, cD — galaxies: nuclei — galaxies: structure 


1. Introduction 

The collisional construction of galaxies from the 
merger of lesser galaxies is thought to be a com¬ 
mon occurrence in the Universe. Coupled with 
the presence of a supermassive black hole (SMBH) 
at the heart of most galaxies (Kormendy & Rich- 
stone 1995; Magorrian et al. 1998, Richstone et al. 
1998), dissipationless mergers have been proposed 
to explain the damaged nuclei in giant elliptical 
galaxies (e.g., Lauer et al. 1995, Faber et al. 1997, 
Rest et al. 2001). Although some galaxy “core¬ 
depletion” is due to the SMBH(s) dining on stars 
that venture to close (e.g., Magorrian & Tremaine 
1999; Zhao, Haehnelt, & Rees, 2002; Yu 2003), 
it is primarily from the gravitational slingshot ef¬ 
fect that the coalescing SMBHs — from the pre¬ 
merged galaxies — have on stars while they them¬ 
selves sink to the bottom of the potential well of 
the newly wed galaxy (Begelman, Blandford, & 


Rees 1980; Ebisuzaki, Makino, & Okumura 1991; 
Makino & Ebisuzaki 1996; Quinlan 1996; Quinlan 
& Hernquist 1997). 

Theory predicts that the orbital decay of two 
such SMBHs should eject a core mass roughly 
equal to the combined black hole masses (Ebisuzaki, 
Makino, & Okumura 1991; Milosavljevic & Mer¬ 
ritt 2001). Current measurements of the central 
stellar deficit are an order of magnitude larger 
than the central SMBH mass, suggesting that 
most elliptical galaxies have undergone multiple 
(«8-10) major-mergers (Milosavljevic & Merritt 
2001; Milosavljevic et al. 2002; Ravindranath, Ho, 
& Filippenko 2002). This result, however, is at 
odds with popular models of hierarchial struc¬ 
ture formation, which predict an average of only 1 
(dissipationless) major-merger event for luminous 
elliptical galaxies (Haehnelt & Kauffmann 2002; 
Volonteri, Madau, & Haardt 2003). 
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Recent advances in our understanding of galaxy 
structure have provided a new framework in which 
to think about, and measure, central mass deficits. 
It is now known that the so-called “power-law” 
galaxies — understood not to have partially de¬ 
pleted cores, nor experienced a (major) dissipa¬ 
tionless merger event — have an undisturbed 
Sersic profile over their entire radial extent. 
That is, their inner light-profiles show no devia¬ 
tion relative to their outer profile (Trujillo et 
al. 2004). On the other hand, the more luminous 
{Mb < —20.5 mag) “core” galaxies display a clear 
flattening of their inner light-profile relative to the 
inward extrapolation of their outer Sersic profile 
(Graham et al. 2003b, 2004; Graham & Guzman 
2003; Trujillo et al. 2004). The ‘break’ in the 
profiles where this transition occurs marking the 
boundary of their relatively unpopulated cores. 

With this new knowledge of what galaxy pro¬ 
files look like, this paper considers the form of the 
core-less galaxy light-profiles when representing 
the original stellar distribution of the disturbed 
profiles, in order to compute the central stellar 
deficit and hence the level of damage to a galaxy’s 
core. 

2. Method & Results 
2.1. Central stellar deficits 

We proceed by quantifying the central stel¬ 
lar deficit as the flux deficit relative to the in¬ 
ward extrapolation of the smoothly curving, stel¬ 
lar distribution outside of any possible, partially 
depleted core. This approach, therefore, does 
not assign any mass deficit to “power-law” galax¬ 
ies, whose inner light-profiles display no clear 
downward deviation from their outer Sersic light- 
profiles. Such a qualitative description can be 
placed on a quantitative footing through employ¬ 
ment of the “core-Sersic” light-profile model (Gra¬ 
ham et al. 2003b), applied in Figure 1 to the “core” 
galaxy NGC 3348. This model consists of an in¬ 
ner power-law and an outer Sersic function. In 
practice (Trujillo et al. 2004), the transition at the 
“break radius” is sharp, providing a 5-parameter 
function^ capable of describing the entire radial 
extent of galaxies with cores. 

^The functional form of the complete core-Sersic model can 
be seen in Graham et al. 2003b and Trujillo et al. 2004. 


The central flux deficit is obtained by differ¬ 
encing the luminosity L{r) = 27rrl(r)dr within 
r = rt, of a) the inwardly-extrapolated outer Sersic 
profile 

Lser(r = rb) = Ierl2TTn ^ Ti{2n, 6 „(r{,/re)^/"), 

ipn ) 

( 1 ) 

and b) a power-law light-profile with slope 7 
matching the observed inner profile slope, and 
intercepting the Sersic model at r = rb, 

I^p-iaw(r' = n) = 2'Krllbl{2 - 7 ). ( 2 ) 

The term Ig is the intensity of the Sersic profile 
at the effective radius Vg enclosing half the galaxy 
light. The exponent ‘1/n’ describes the curvature 
of the light-profile, and is simply a function of 
n such that r( 2 n) = 2 ri( 2 n, 6 „), where F and Fj 
are the complete and incomplete gamma functions 
respectively. The intensity at rb is denoted by /&. 

This procedure has been applied here to the 7 
bona fide core galaxies from Trujillo et al. (2004)^. 
The resultant flux deficits are 0.07-0.7% of the 
total galaxy flux. Allowing for the fact that 
these galaxies are at a range of distances from us, 
the apparent magnitude differences m{r = rb) = 
—2.51og[Lser — 7vp_iaw] Were converted into abso¬ 
lute magnitudes (M) using the galaxy distance es¬ 
timates from Tonry et al. (2001) and then con¬ 
verted into units of solar flux using an absolute 
i?-band magnitude for our Sun of Mq =4.46 mag 
(Cox 2000). Finally, these values were trans¬ 
formed into solar masses assuming a stellar mass- 
to-light ratio of 3.0 (Worthey 1994). Such a ra¬ 
tio is representative of an evolved (i.e., faded) 
12 billion year old, single stellar population ob¬ 
served through an i?-band optical filter. The mass 
deficits Mdef are given in Table 1. 

2.2. Constraints on the number of dry 
mergers 

Models of hierarchical structure formation pre¬ 
dict that galaxies will collide; indeed, this phe¬ 
nomenon has been observed for many years. The¬ 
oretical expectations for the ejected core mass, af¬ 
ter the violent unification of galaxies containing 

^Clarifying note: The structural parameters tabulated in 
Trujillo et al. (2004) were obtained using the approximation 
bn ~ 1.999n—0.3271, and the values are for each galaxy’s 
outer profile as if it had no partially depleted core. 
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SMBHs, scale as 0.5-2A^Mbh, where Mbh is the fi¬ 
nal BH mass and N is the number of merger events 
(Milosavljevic et al. 2002). The variable term at 
the front is because equal-mass mergers scour out 
more stars than a collision involving a lesser-mass, 
secondary galaxy. Thus, by knowing a galaxy’s 
central stellar deficit (Table 1), and its black hole 
mass Mbh, one can place constraints on the extent 
of this merger process. 

We have estimated the SMBH mass of each 
galaxy using two techniques. First, we employed 
the Mbh-u relation of Gebhardt et al. (2000), in 
which Mbh is derived from the galaxy velocity dis¬ 
persion (T. Using the steeper Mbh-u relation of 
Ferrarese & Merritt (2000) had little difference be¬ 
cause of the different zero-points in these two re¬ 
lations. The second SMBH mass estimate came 
from the independent Mbh“U relation (Graham 
et al. 2001, 2003a). This relation is as strong 
as the Mbh-u relation and has the same small 
degree of scatter, but has the advantage that it 
is a purely photometric technique, requiring only 
(uncalibrated) images rather than (telescope-time 
expensive) galaxy spectra. Figure 2 shows each 
galaxy’s depleted core mass plotted against both 
estimates of the SMBH mass. The ratio of the 
stellar mass deficit to central black hole mass has 
a mean ± average deviation of 2.4 ± 0.7 (panel a, 
Mbh-u) and 2.1 ±1.1 (panel b, Mbh-u), consistent 
with these galaxies having experienced one major 
(i.e., equal-mass) dry merger event. 

3. Discussion 

Support for the above mass deficits comes from 
the agreement with the merger simulations of 
Makino & Ebisuzaki (1996), the (cusp regener¬ 
ation) hierarchial merger models of Volonteri et 
al. (2003), and the theoretical expectations of 
Haehnelt & Kauffmann (2002). Although the lat¬ 
ter Authors predicted a median number of equal 
mass mergers for faint (power-law) and bright 
(core) elliptical galaxies of 1 and 3 respectively, the 
number of major mergers since the last collision 
that involved substantial gas accretion (Mgas > 
Mbh) is 0 and 1 respectively. The presence of gas is 
important because it dilutes the wrecking ball ac¬ 
tion of the SMBHs on the stars because it, rather 
than the stars, fosters the coalesence of the black 
holes and it can lead to the creation of new stars 


(see, e.g., Zhao et al. 2002, and references therein). 
Therefore, our conclusions that power-law galaxies 
do not have partially depleted cores from galaxy 
collisions, and that the number of (dissipationless) 
major mergers producing luminous galaxies, with 
Mr ~ —22.5 mag, is equal to about 1, are sup¬ 
ported by current cold dark matter models of hi¬ 
erarchial structure formation. 

Given Mdef ~ Mbh, one may wonder if some 
depleted cores might have formed from the run¬ 
away merging of stars (e.g., Begelman & Rees 
1978; Quinlan & Shapiro 1989) within what may 
once have been the dense cores of massive ellipti¬ 
cals (Graham & Guzman 2003), rather than from 
the scattering of stars from coalescing black holes. 
The first objection to such a process would be 
those cases in which Mdef is actually greater than 
Mbh- Such a mechanism would also require a cer¬ 
tain level of refinement, such as re-populating the 
loss-cone, in order to explain the absence of (re¬ 
solved) cores in less luminous elliptical galaxies. 
The expected break radii — derived from a cen¬ 
trally depleted Sersic model — are not observed 
amongst the “power-law” galaxies. For example, 
if the 12 “power-law” galaxies in Trujillo et al. 
(2004) had cores with inner power-law slopes rang¬ 
ing from 0.0-0.3, then, assuming Mdef = Mbh, 
they should have break radii of 0.17-0.5 arcsec- 
onds, which they do not. 

It is pertinent to inquire why previous esti¬ 
mates of the central mass deficit are larger than 
the values obtained here. One reason is that it 
had been assumed that every galaxy once had a 
steep isothermal p{r) ~ r~^ core before any merg¬ 
ing black holes wreaked their havoc. There is, 
however, no observational evidence that any such 
universal density-profile p(r) exists, or once ex¬ 
isted amongst the “power-law” galaxies. In actn- 
ality, a luminosity-dependent range of inner pro¬ 
files shapes is now known to exist (Gebhardt et al. 
1996; Graham & Guzman 2003; Balcells, Graham, 
& Peletier 2004). No single power-law slope can 
be used to approximate the initial, undisturbed, 
stellar distribution of all elliptical galaxies. 

Milosavljevic et al. (2002), for example, had 
defined the onset of partially depleted cores as 
the radius where the negative, logarithmic gradi¬ 
ent of the deprojected light-profiles (i.e., the spa¬ 
tial, luminosity-density profiles) equaled 2. That 
is, where the observed slope matched that of the 
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isothermal model. They then derived central, stel¬ 
lar mass deficits from the difference between the 
observed density-profile and the inward extrapola¬ 
tion of the isothermal model from this point. This 
approach to estimating the depleted core mass is 
illustrated in Figure 3, where we show both the ob¬ 
served light-profile and the (deprojected) density- 
profile of the core-less, “power-law” galaxy NGC 
5831. Such a prescription encounters a number of 
difficulties. An extraordinary level of fine-tuning 
would be required to deplete stars over the full ra¬ 
dial extent of the initial isothermal core, assuming 
one existed, but not beyond the final core radius, 
which one assumes still has its original slope of -2 
today. Moreover, with such an approach, many 
of the “power-law” galaxies which are not (tradi¬ 
tionally) recognized as having a partially depleted 
core will be assigned one. Such theoretical core 
radii (and mass deficits) are not only questionable, 
but excessively large — sometimes greater than 1 
kpc — and don’t match the observed break radii 
in core galaxies, which are invariably less than a 
few hundred parsecs (see Table 1). 

A second approach to estimating the central 
deficit had been to use break radii derived from the 
Nuker model (Lauer et al. 1995), and to assume 
an isothermal core once existed inside of this ra¬ 
dius. With such an approach the estimated central 
deficits are again too high because of a) the exces¬ 
sively steep isothermal model that is assumed, and 
b) because, as described in Graham et al. (2003b) 
and shown in Trujillo et al. (2004), Nuker-derived 
break radii typically overestimate the actual break 
radii by factors of 2 to 5. 

The mass deficits obtained with the above two 
methods are 3-30 times larger than our values, 
and are shown in Figure 2 for comparison. Our 
new measurements of the mass deficit are signifi¬ 
cantly (99.4%, from both a K-S test and Students’ 
T-test) different to those values obtained using the 
isothermal assumption, and reveal that the galac¬ 
tic merger history of the Universe, at least for 
massive elliptical galaxies, is roughly an order of 
magnitude less violent than previous observational 
analyses (Milosavljevic & Merritt 2001; Milosavl- 
jevic et al. 2002; Ravindranath, Ho, & Filippenko 
2002 ) had suggested. 

We plan on analyzing more galaxies, and mea¬ 
suring the ellipticity of their evacuated core re¬ 
gion. This may shed light on the orientation of 


the initial orbits of the black holes, and allow one 
to explore any possible correlation with the host 
galaxy ellipticity. A greater range of data will also 
allow one to explore whether there is any trend 
between Mdef/Mbh and galaxy magnitude. Such a 
correlation may be expected if bigger “core” galax¬ 
ies have experienced more dissipationless merger 
events than less luminous “core” galaxies. 

Support for proposal HST-AR-09927.01-A was 
provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by 
the Association of Universities for Research in As¬ 
tronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 1.— Observed, major-axis, i?-band, sur¬ 
face brightness proHle of NGG 3348. The solid 
line is the best-fitting core-Sersic model, while the 
dashed line is the best-Htting Sersic model to the 
data beyond the break radius rt,=0.35 arcseconds. 
The flux deficit is illustrated by the area desig¬ 
nated as the ‘depleted zone’, corresponding to a 
mass deficit of 300 million solar masses. Data 
points from Trujillo et al. (2004), and supercede 
those shown in Graham et al. (2003b). 



Fig. 2.— Central mass dehcits evaluated from 
three techniques versus the central black hole 
masses derived using a) the galaxy Sersic index 
n (Graham et al. 2001, 2003a), and b) the veloc¬ 
ity dispersion a (Gebhardt et al. 2000). 6-pointed 
stars: from Milosavljevic et al. (2002) using the 
method illustrated in Figure 3b. 5-pointed stars: 
derived using Nuker-model break-radii and equa¬ 
tion 41 from Milosavljevic & Merritt (2001). Filled 
circles: derived here using the logic illustrated in 
Figure 1. The solid line shows a one-to-one rela¬ 
tion, the dashed line shows Mdef = 2Mbh. Typical 
errors for the points marked with a filled circle are 
roughly log 2 along both axes, stemming from a 
~ 20% and ~ 15% uncertainty on the value of 
n and a respectively, and from proHle Htting and 
galaxy distance uncertainties. 
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Table 1 
Galaxy data 


Galaxy 

a 

(km s“^) 

Mgal 

(i?-mag) 

n 

core-Sersic 

< (pc) 

Mdei 

Isothermal 
r/ (pc) Mdef 

Nuker 

r/ (pc) Mdef 

NGC 2986 

268 

-22.49 

5.28 

0.69 (97) 

7.02 

2.9 (410) 

35.5 

1.24 (174) 

26.7 

NGC 3348 

237 

-22.76 

3.81 

0.35 (70) 

3.01 

2.5 (490) 

66.1 

0.99 (198) 

26.5 

NGC 4168 

186 

-21.92 

3.12 

0.72 (108) 

1.21 

8.3 (1250) 

38.0 

2.02 (303) 

23.5 

NGC 4291 

284 

-21.61 

5.44 

0.37 (47) 

5.57 

1.3 (170) 

17.8 

0.60 (076) 

16.7 

NGC 5557 

253 

-23.08 

4.37 

0.23 (51) 

2.13 

2.0 (440) 

41.7 

1.21 (269) 

39.6 

NGC 5903 

210 

-22.69 

5.09 

0.86 (141) 

8.44 

4.7 (780) 

38.0 

1.59 (262) 

25.5 

NGC 5982 

250 

-22.99 

4.06 

0.28 (57) 

3.12 

2.2 (450) 

39.8 

0.74 (151) 

23.5 


Note. —Column 1: New General Catalog (NGC) numbers. Column 2: Velocity dispersions a from Hyper¬ 
cat (littp://www-obs.univ-lyonl.fr/hypercat/). Column 3: Absolute i?-band galaxy magnitudes Mgai derived 
from the best-fitting (sharp transition) core-Sersic model; the parameters from which, including the Sersic 
index n and the break radii ri, (in arcseconds and parsecs), can be found in Trujillo et al. 2004. The asso¬ 
ciated central mass deficits Mdef (in units of 10 ® solar masses) have been derived from the difference in flux 
between the observed light-prohle and the inward extrapolation of the outer Sersic profile. The ‘Isothermal’ 
quantities have come directly from Milosavljevic et al. 2002 who used the technique illustrated in Figure 3b. 
The ‘Nuker’ break radii have come from Rest et al. 2001, and the Nuker-derived mass deficits were obtained 
using equation 41 from Milosavljevic & Merritt 2001 with the required 7 values taken from Ravindranath et 
al. 2002 . 
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Fig. 3.— a) Observed, major-axis, i?-band, sur¬ 
face brightness profile of the core-less galaxy 
NGC 5831 (data taken from Trujillo et al. 2004). 
A 3-parameter Sersic model (solid line) adequately 
describes the stellar distribution on both the nu¬ 
clear (< 1") and global scale; the dashed line 
shows the model extrapolated beyond the data. 
No characteristic downward break in the inner 
light-profile, that would signify damage caused 
by merging SMBHs, is evident, b) Spatial den¬ 
sity profile (curved line) of NGC 5831 obtained 
by deprojecting the light-profile model in panel 
a). Assuming an isothermal model, p{r) ~ 
once existed (straight line), one can assign a the¬ 
oretical core radius Vf, — where the logarithmic 
density profile has a slope of -2 — and a central 
mass deficit. Milosavljevic et al. (2002) obtained 
rb=220 parsecs and M^ef = 5.6 x 10® Mq. 
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